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Abstract. We consider here square tilings of the plane. By ex- 
tending the formalism introduced in [3] we build a correspondence 
between plane maps endowed with an harmonic vector and square 
tilings satisfying a condition of regularity. In the case of periodic 
plane square tiling the relevant space of harmonic vectors is actu- 
ally isomorphic to the first homology group of a torus. So, periodic 
plane square tilings are described by two parameters and the set 
of parameters is split into angular sectors. 

The correspondence between symmetry of the square tiling and 
symmetry of the plane maps and harmonic vectors is discussed and 
a method for enumerating the regular periodic plane square tilings 
having r orbits of squares is outlined. 



1. Introduction 

The face-to-face plane square tiling [0, 1]^ -|- is the most basic 
example of plane tiling. An interesting variation, found on many city 
squares, is the plane tiling by squares of two different size and no face- 
to-face adjacency (see Figure [T]). The next interesting case, i.e. the 
classification of plane tilings with no face-to-face adjacency and 3 orbits 
of squares of different square sizes was begun in [7] and completed in 
[HI Uni H] by using classical methods of tiling theory. In parallel to this 
in [3] an electrical network formalism was introduced for the dissection 
of rectangles into squares of unequal length (see |6j for a very good 
history of the problem). 

In Section [2| we propose a combinatorial formalism for describing 
plane square tilings satisfying the condition that there is no line con- 
taining an infinity of edges of square. The correspondence associates 
to any such tiling a plane map Ai and an harmonic vector w on it. 
Reversely, for every plane map and harmonic vector we define an un- 
extendible plane square packing, which is actually a tiling in many 
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Figure 1 . The city square tiling and the corresponding 
plane map 



cases. This directly extends the work of [3] for dissection of rectangles 
to square plane tilings. 

In Section |3| we identify, in the periodic case, the harmonic vectors 
to homology classes. That is, to any edge finite periodic tiling T of 
the plane by squares, we associate a periodic plane maps Ai. Then the 
tiling T is described in terms of two parameters which correspond to 
the first homology group of Ai quotiented by its group of translations. 
However, in the non-periodic case, such an identification is not possible 
and the variety of tilings is much greater. 

In Section |4| we see that the space of possible parameters of a square 
packing does not correspond in general to the full homology group. In 
general, the two dimensional space of the first homology group of a 
torus is split into angular sector containing the origin. We look at how 
the 8 periodic square tilings of [U [101 [1] cire thus described. 

The wallpaper symmetry groups of such tilings are restricted and in 
Section [5] we consider how the symmetries of a periodic square tiling 
correspond to the symmetries of the corresponding toroidal map. We 
then introduce a systematic method for enumerating the periodic reg- 
ular square tilings with r orbits of squares. 

2. Harmonic vector formalism 

For a tiling of the plane or a rectangle by squares aligned along the 
X axis, we call horizontal edge, respectively vertical edge, the edges of 
squares aligned with x-axis, respectively y-axis. An horizontal face 
is the set of horizontal edges contained in an horizontal line. See an 
example on Figure [2j Similarly we define vertical faces. 
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Figure 2. An horizontal face 
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Figure 3. Dissection of rectangle into squares and cor- 
responding electrical network 

In [3j the problem of dissecting a rectangle into squares of unequal 
size was considered. There, to any such dissection into n squares is 
associated an electrical circuit on n wires each of conductance 1. This 
electrical circuit forms a finite planar graph between between two points 
A and B. For an horizontal face, the sum of lengths of above edges 
should be equal to the sum of lengths of below edges; this is translated 
in conservation of currents at vertices. Similarly equality of sum of 
lengths at vertical edges correspond to voltage equality in faces. Picture 
[sjfrom [3] indicates the correspondence in an example. 

For tiling of the plane by squares, the correspondence goes along the 
same lines but some technical conditions are needed. A square tiling is 
called edge finite if all its vertical and horizontal faces contain a finite 
number of vertical and horizontal edges. A vertex is called singular if 
it is contained in exactly 4 different squares. A plane square tiling is 
called regular if it is edge finite and has no singular vertices. 

By a plane graph Ai we mean a graph which can be drawn on the 
plane such that any two edges if not disjoint interest only at their 
vertices. The faces of Ai are the connected components of the plane 
minus vertices and edges. We denote by V{J^), E{Jli) and F(A^) the 
set of vertices, edges and faces of A^. A plane map is a connected plane 
graph such that any face is bounded and any compact subset of the 
plane contains only a finite number of vertices and edges. Such a map 
is necessarily infinite. A plane map is called periodic if it is invariant 
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under translation in two non-parallel directions. Of course periodic 
plane maps correspond to toroidal maps and we can use the notions of 
vertices, edges and faces for tori. 

Given a plane map Ai a directed edge "e^ is a pair (f , e) with v a 
vertex and e an edge containing v. For a vertex v, edge e, or face / we 
denote by ^(f), ^(e), or ^(/) the set of directed edges originating 
from V, contained in e, or having / on the right. An harmonic vector 
w = (w-^) of a plane map M. is a vector defined over the set of directed 
edges such that for any v G V{M), e G E{M), f G F{M) we have: 

= = ^ = w^. 

-te'^iv) "^e^{e) "^e^(/) 

We denote by Harm{Ai) the vector space of harmonic vectors of A^. If 
is a periodic map then we denote by HarmP^'^{A4.) the vector space 
of periodic harmonic vector. For any edge e = (1^1,^2) we have two 
directed edges (wi,e) and (f2,e) and the relation W(^^^^e) + w;(^2,e) = 0. 
Thus we can define = |ti'(t,j,e)| = |iy(t,2_e)|- When drawing plane maps 
and their harmonic vectors, it is convenient to choose an orientation 
on each edge and to write down the value of the harmonic vector next 
to it. 

We now construct to an edge finite square tiling T a plane map Ai 
and a vector w G Harm{A4). In this correspondence, every square 
corresponds to an edge. Faces of Ai correspond to horizontal faces of 
T and vertices of Ai to vertical faces of T. For regular square tilings 
the map Ai is uniquely defined. However, if T has a singular vertex 
then there are two ways to define the horizontal and vertical faces and 
this affects the map Ai: either the singular vertex is contained in one 
horizontal face and two vertical faces from above and below or it is 
contained in one vertical face and two horizontal faces from left and 
right (See in Figure |4] an example of this ambiguity). Note that all the 
plane tilings by squares considered in [7j are regular. 

Theorem 1. For an edge finite plane square tiling T , there exists a 
plane map Ai and a harmonic vector w G Harm{Ai) such that: 

(i) Every square ofT corresponds to an edge e of Ai with being 
the size of the square. 

(ii) For a choice of horizontal and vertical faces in T the map Ai 
is uniquely defined; its vertices correspond to horizontal faces of 
T and its faces to vertical faces of T ■ 

(iii) // T is periodic and the choice of horizontal and vertical faces 
is done in a periodic way then the map Ai and the vector w are 
also periodic. 
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Figure 4. A singular vertex and the two possible graphs 
Ai locally around it according to the choice of vertical 
and horizontal faces 

(iv) // T is regular then Ai is determined uniquely. In that case, if 
T is rotated by 90 degrees then M. is changed into its dual Ad* . 

Proof. Let us take a choice of horizontal and vertical faces in T. We 
put a vertex in the middle of every horizontal face. Also, for every 
square of T we put an edge between the upper horizontal face and the 
lower horizontal face, see Figure |4] for an example. Clearly, the vertical 
faces of T correspond to faces of and (ii) follows. If C is a square 
of T then there are vertices f , v' corresponding to the horizontal faces 
in which the upper, respectively lower, edge of C are contained. We 
orient ~^ from v to v' and set w-^ = d with d the size of C. For the 
reversed directed edge, we take the opposite value. One checks that 
the vector w thus defined is harmonic and (i) follows, (iii) is clear since 
the above process is uniquely defined once the horizontal and vertical 
faces are chosen. If T is regular then there is only one possible choice 
of horizontal and vertical faces. Furthermore, in that case rotating by 
90 degrees corresponds to taking the dual and (iv) follows. □ 
We now want to define a square tiling of the plane from a plane map 
and an harmonic vector. Before that we state a lemma that is essential 
to this construction: 

Lemma 1. If M. is a plane map and w G Harm{Ai) then: 

(i) If V is a vertex of M. and ei, . . . , are the directed edges 
coming from v cyclically ordered then up to rotation the sign 
shape is (+", 0", — 0"^) with a + n + b + m = N 

(ii) If f is a face of M. and ci, . . . , cn are the directed edges 
around f cyclically ordered then up to rotation the sign shape 
is 0", -^ 0"") with a + n + b + m = N 
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Proof. This property is established in Lemma 4.11 of [3] in the context 
of dissection of rectangles by squares. One notices that the proof given 
there does not depend on the hypothesis of conductance equal to 1. 

Now for a given plane map Ai and w G Harm{A4), let us select a 
vertex v of A4. We can find a block B of faces that contains v and cut 
all edges outside of B. Now by adding edges with the right conductance 
around B one can close the set of circuit and get a set of wires with the 
current coming from a vertex A and arriving at B. Thus we can apply 
Lemma 4.11 in [3J and conclude. □ 

We can encode dissection of rectangle by rectangles by a similar 
formalism as the one of [3J. In the correspondence, rectangles of side 
length a, b are associated to wires of conductance a/b. Thus the spirit 
of the above lemma is that a finite set of squares can be completed to 
a rectangle by adding some rectangles. 

For a plane map Ai the medial map Med{Ai) is the map obtained by 
putting a vertex on every edge of M and having two edges adjacent if 
and only if they share a face and a vertex. The resulting map Med{Ai) 
is 4-valent and we have Med{M) = Med{M*). The dual {Med{M))* 
is bipartite and its bipartite components correspond to and M* . 

For a general map M. and w G Harm{M.) we cannot guarantee that 
we will get square tiling of the plane. Let us call an unextendible square 
packing T, a plane square packing such that around every square there 
are only squares and no free space. An example of unextendible square 
packing which is not a tiling is obtained by taking the subdivision of 
the upper half plane M+ x M into unit squares. Then one subdivides 
the squares in [0, 1] x M into squares of sides 1/2 and then the squares 
in [0, 1/2] X M into squares of sides 1/4. Continuing in this process one 
gets an unextendible square packing which is not a tiling. 

Theorem 2. (i) If M. is a plane map and w G Harm{M.), w ^ Q 
then there exists an unextendible square packing T{Ai,w) such that by 
applying the construction of Theorem [I] to the corresponding choice of 
horizontal and vertical faces we get back Ai and w. 

(a) The unextendible square packing T{Ai,w) is unique up to trans- 
lations and is simply connected. 

Proof. By applying Lemma [T] we get that the signs around a vertex or 
a face of Ai are organized in a block of + followed by a block of — with 
possibly between them. This gives a simple way to reconstruct the 
square tiling from Ai and w. Simply take one edge e of and assign 
it to a square C of opposite corners (0, 0) and {we, We). Since e belongs 
to two vertices this allows to define the horizontal faces over C and 
below it. The edge e is also contained in two vertical faces and those 
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define vertical edges on the left and right of C. Let us call two edges of 
Ai adjacent if they share a vertex and a face. So, once the position of 
C is chosen the position of the squares corresponding to edges adjacent 
to e is also set. After having put the squares around C we can continue 
in this way for other edges and so we get the uniqueness of T{A4,w) up 
to the position of the initial point. But in order to prove the existence 
we need something else. We need to prove that by defining the position 
iteratively in the graph Med{M.) we always get coherent positions. 

Let us take a cycle C = {ei, . . . , cat = ei} of edges with Cj adjacent 
to Cj+i. We associate a square Ci to every edge Cj according to the 
above defined procedure and we need to prove that Ci = Cn- Without 
loss of generality we can assume that C does not self-intersect. Let 
us now denote by Fi, . . . , Fm the faces inside of C. By applying 
the relations satisfied by w around vertices, we can modify C so that 
around any vertex of C the edges turn inward of it. If we change the 
cycle C by removing one face Fi which touches the boundary then we 
do not change the position of the square. This is because the vector 
w is harmonic and the sum of edge weights is zero around every face. 
Thus we can remove faces Fi one by one until one is reduced with the 
trivial cycle {ei} which proves that Cn = Ci. Assertion (i) and (ii) 
follow. □ 



3. Homology groups 

In this section we consider the space of harmonic form, first in the 
case of a periodic plane graph: 

Theorem 3. (i) For a periodic plane map Ai and w G Harm^'^^{M.) 
the unextendihle square packing T{A4,w) is a periodic square tiling of 
the plane. 

(ii) For a periodic plane graph Ai, the space HarmP^'^ {J^) is of 
dimension 2 and is isomorphic to the first homology group of a 2- 
dimensional torus. 

Proof. Assertion (i) follows from the uniqueness condition of Theorem 
|2} The fact that it is a tiling follows from the fact that 1^(6)1 > C 
with C > or |w(e)| = 0. This prevents the formation of accumulation 
points and guarantees that T{M.,w) is a tiling. 

Denote by L = TjVi + Zt'2 the group of translation of Ai with Vi not 
collinear to V2. The torus T = A^/L is a map with vertices, edges and 
faces, i.e. it is a cell complex. Moreover, the space HarmP^'^ {Ai) is 
isomorphic to Harm{T). 
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For an edge e = (f i, ^2), the condition ^-^g;^^^^ w-^ = is rewritten 
as W(t,i,e) + wi^v2,e) = 0. This simply means that if we choose an orienta- 
tion on every edge then the space HarmiT) is defined over M'^C'''). We 
can thus build a chain complex: 

If one sets the standard scalar products on ]R^*^''\ M.^^'^^ and M.^^'^^ 
then we can define the adjoint and (ij- The space HarmiT) is then 
Ker((i2) H 'Kei{d\) and by Hodge Theory (See, for example, Theorem 
7.55 in [9]) it is identified with the first homology group ifi(T, M) with 
coefficients in M. Since T is a torus, ifi(T,]R) is of dimension 2. □ 
As a consequence of (ii) we know that if Ai is periodic along a lattice 
L and w is an harmonic vector periodic along a lattice L' G L then 
w is periodic along L as well. If Ai is periodic then let us write the 
lattice L of translations asL = Zti + Zt2- The translations along 
t i correspond to a cycle 7i in the quotient A4/L. Let us write 7, as a 
cycle of vertices 

{vlvl...,v[} with vl = vl + ~t, 

and the pair = (ff , ff^^) being edges. We can now define the integral 
of w G HarmP^'^{M.) as 

i-i 



fc=l 

The square tiling T(A^,w) is periodic and we can write its lattice of 
translation as L' = Z t [ + Z t 2- The y component of t is w. For 
the X component, one must proceed differently. The harmonic space 
HarmP^^{Ai*) is isomorphic to HarmP'^'^{M.) and one can realize the 
isomorphism explicitly. It suffices to turn the oriented edges of M. by 
90 degrees and assign to them the same edge weight as in If w* 
is the induced harmonic vector and 7* are the induced cycles then the 
X component of t ^ is equal to w* . See in Figure j5j a complete 
example. 

Theorem [s] (ii) is classical and corresponds to the identification of 
harmonic forms on a compact manifold with real homology classes. 
However, it is well known that Hodge Theory collapses in the non- 
compact case. This allows us to build many interesting square tilings: 

Theorem 4. There exists a square tiling of the plane with tiles of 
size all different such that the size of the tiles grows linearly with the 
distance from the origin. 
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Figure 5. A periodic square tiling with translation vec- 
tors t[ = {b — a,a + b + c) and t'2 = (—a — b,a — b). See 
also the associated plane map Ai and its dual Ai* 

Proof. Let us consider the plane map formed by a family of squares. 
Let us assign a value a{k,l), b{k,l) to horizontal, respectively vertical 
edges according to Figure |6j The fact that this defines an harmonic 
vector corresponds to the following equations for fc, / G Z: 

( a{k,l)-a{k-l,l) + b{kj)-b{k,l-l) = 0, 
\ a{k,l) + b{k + l,l)-a{k,l + l)-b{k,l) = 0. 

The two solutions below 

wi = {a{k, /) = 1 and b{k, I) = 0}, 
W2 = {a{k,l) = and b{k,l) = 1} 

form a basis of HarmP'^'^{7?). But Harm{7?) is not reduced to it. The 
next most interesting harmonic vectors are: 

= {a{k, I) = I and b{k, I) = k} 
W4 = {a{k,l) = k and b{k,l) = — /} 

Now if we set the harmonic vector 

w = (l/4)wi + (1/3)^2 + W3 + V2wi 

then one sees that \w{e)\ = \w{e')\ implies e = e' and that the growth 
of |w(e)| is linear with the distance from the origin. Since We goes to 
infinity on each horizontal and vertical direction of 7?, we get that the 
unextendible square packing T(Z^,w) is actually a tiling. □ 
The above theorem could be seen as a counterexample to a problem 
posed by C. Pomerance and mentioned in [7]: does there exist a square 
tiling of the plane with all squares of different size such that after 
ordering the sizes the growth is less than exponential? However the 
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Figure 6. The edge weight of the map 1? 

question is not exphcitly stated there and if one imposes that square 
sizes are integral then we do not have the answer to the question. 

4. Types of square tilings 

Given a periodic plane map and w G Harm^'^^'{Ai) we have a 
corresponding periodic square plane tiling l~{Ai,w). By Theorem Is] 
we can find a basis {w^,w'^) of HarmP'^^{AA). So, w = aiw'' + a2Vu 
and this allows us to define the corresponding 5'g-domains: 

Definition 1. Given a periodic plane map Ai and a basis {w^jw"^) of 
Harm^^^ (Ai) we define: 

(i) Plane{M. , {w^ , w'^}) to he the plane of the parameters (0:1,02). 

(ii) For an edge e of M., we choose an oriented edge "e^ G ^(e) and 
define a line L(e) to he the set of (Q;i,a;2) £ Plane{Ai, {w^,w'^}) such 
that aiw^ + a2'W^ = 0. 

(Hi) A Sq-domain to he a connected set D of Plane{M, {w^jw"^}) 
which is delimited hy lines L(e) such that aiw^ + a^w"^ is of specified 
sign e(e) for a = (ai, 0^2) G D. 

Given Ai and [w^jw"^) it is very easy to determine all possible Sq- 
domains. For every edge e of we associate the line aiW^+a2W?^ = 
in the plane (ai, 02) • All such lines pass by the origin and so the plane 
is partitioned in a number of angular sectors which are exactly the 
5'g-domains. 

Take now a point near a line L{e). If one goes towards L{e) then 
the length w{e) of the corresponding square vanish. After one passes 
through L{e), the square reappears but in a different position. In Figure 
[7] the eight unilateral plane tilings with three sizes of faces are listed. 
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They are all regular and we give their corresponding plane maps Ai. 
One notices that the square tilings (2), (3) and (4) are associated to 
the same periodic plane map. 

5. Symmetry and classification considerations 

Suppose that we have a periodic plane map with a group T. Any 
/ G r defines an action on H anmP^"^ {M) . The induced action of F 
is actually the action of the point group of V. 

For a square tiling, rotation axis of order 3 or 6 are forbidden and 
this leaves pi, p2 and p4 as possible rotation groups. However, we do 
not know which of the 17 wallpaper groups can occur as group of a 
periodic square tiling. 

Let us now see how groups of M. and their square tiling T(A^,w) 
are related: 

Proposition 1. (i) If f is a symmetry of a periodic plane map Ai and 
w G Harm^^^{Ai) the f induces a symmetry of T{Ai,w) if and only 
if f*iw) = ±w. 

(a) If has symmetry p2 then all square tilings T{M.,w) also have 
symmetry p2 at least. 

(Hi) Ifw is chosen at random in Plane{Ai, {10^,10"^}) then T{Ai,w) 
has symmetry pi, p2 or p4. 

Proof, (i) is quite clear. The symmetry p2 is the antipodality so of 
course it satisfies condition (i) and maps itself to the tiling. For a 
random w G Harm^^'^{M.) no refiection can preserve w and so the 
symmetry of a generic square tiling is pi, p2 or pA. □ 

In order to explain the symmetry p4 one has to introduce the notion 
of self-duality. If one looks at the definition of the harmonic space 
Harm{A4.) one remarks that vertices and faces play identical roles. 
A map is called self-dual if the dual Ai* is isomorphic to Ai. For 
the square tiling T{Ai,w) taking the dual corresponds to rotating by 
90 degrees. Thus if a regular square tiling has a 4-fold axis then its 
corresponding plane map Ai is necessarily self-dual. 

Figure [Tjgives for each of the eight plane tiling with 3 orbits of squares 
of different size the corresponding plane map AA.. The first five such 
tilings have symmetry p2 and are defined by two parameters. The last 
three also have a line of symmetry and this comes from a symmetry of 
the plane map and harmonic vector. If a symmetry preserves a vector 
then this fixes some parameter and thus those tilings are defined by 
just 1 parameter. 

We cannot use the construction that we have introduced for the 
enumeration of plane square tilings with r orbits of squares. The first 
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(3): p2 (4): p2 




(7): pgg (8): pgg 



Figure 7. The 8 square plane tiling with 3 orbits of 
squares and their corresponding plane maps 
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reason is the restriction of edge finiteness. The second reason is that if 
there is a singular vertex then there is an element of choice in the con- 
struction of the map Ai, which greatly complexifies the enumeration. 

However, if one restricts to regular square tilings, then we can de- 
sign a possible enumeration strategy. If the square tiling T does not 
have a 4-fold axis of symmetry but has r orbits of squares then the 
corresponding map A4 has r orbits of edges. But if T has a 4-fold 
axis of symmetry then it comes from the self-duality of A4. The me- 
dial operation satisfies Med{Ai) = Med{Ai*) and any self-duality of 
Ai corresponds to an ordinary symmetry of Med{Ai). Moreover, the 
maps Ai and Ai* can be reconstructed easily from Med{Ai) and this 
was already used for enumeration of self-dual maps in [5]. 

Thus we have to enumerate the toroidal maps having r orbits of 
vertices in order to get the regular square tilings with r orbits of squares 
by applying the technique of Delaney symbol [1] and the methods of 
exhaustive enumeration P]. This should actually be possible for small 
values of r like r = 4, 5. However, as noted in [7J what would be really 
nice is if the existing square plane tilings of Figure [7] were actually used 
in city squares. 



References 

[1] A. Bolcskei, Classification of unilateral and equitransitive tilings by squares of 
three sizes, Beitrage zur Algebra und Geometrie 41 (2000) 267-277. 

[2] G. Brinkmann, Isomorphism rejection in structure generation programs, Dis- 
crete mathematical chemistry, DIMACS Ser. Discrete Math. Theoret. Comput. 
Sci. 51 (Amer. Math. Soc, 2000) 25-38. 

[3] R. L. Brooks, C. A. B. Smith A. H. Stone and W. T. Tutte, The dissection of 
rectangles into squares, Duke Math. J. 7 (1940) 312-340. 

[4] O. Delgado, D. Huson and E. Zamorzaeva, The classification of 2-isohedral 
tilings of the plane, Geometriae Dedicata 42 (1992) 43-117. 

[5] M. Dutour Sikiric, M. Deza, ^-regular and self-dual analogs of fullerenes. Math- 
ematics and Topology of Fullerenes, Lecture notes in physics, 2010, to appear. 

[6] I. Gambini, Quant aux carres carreles, PhD thesis Universite de la 
Mediterranee, 1999. 

[7] B. Griinbaum and G. C. Shephard, Tilings and Patterns, W. H. Freeman 
(1987) 72-81. 

[8] H. Martini, E. Makai and V. Soltan, Unilateral tilings of the plane with squares 
of three sizes, Beitrage zur Algebra und Geometrie 39 (1998) 481-495. 

[9] M. Nakahara, Geometry, Topology and physics, Institue of physics publishing, 
1989. 

[10] D. Schattschneider, Unilateral and Equitransitive tilings by squares. Discrete 
Comput. Geom. 24 (2000) 519-515. 



14 



MATHIEU DUTOUR SIKIRIC 



Mathieu Dutour Sikiric, Rudjer Boskovic Institute, Bijenicka 
10000 Zagreb, Croatia, Fax: +385-1-468-0245 
E-mail address: mdsikir@irb.hr 



